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Abstract 

In this sudy, a dispersion measure for stable distributions is pro- 
posed by measuring the 'average curvature' of the characteristic func- 
tion of a stable random variable. The 'average curvature' is to be 
called the characteristic curvature and it resembles the concept of 
Fisher information. The intellectual motivation comes from the un- 
certainty principle for Fourier transform pairs. Characteristic curva- 
ture reduces to that of standard deviation in the Gaussian case and 
possesses great many analytical features, for example the scaling be- 
haviour of the type t * and its explicit connection to the Euler Gamma 
function. It is also directly related to the Hurst exponent H and the 
fractal dimension 2 — H. 

1 Introduction 

Stable distributions is an interesting and important class of probability dis- 
tributions. They were discovered explicitly by Paul Levy in 1925 pQ. They 
possess many interesting properties, most importantly they are by definition 
invariant under addition, up to a scale. Noteworthly they have power-law 
type of decay and therefore they are an excellent model for modelling many 
natural phenomena, such as earthquakes, financial returns, and a multitude 
of social phenomena such as size distributions of cities and firms [6]. The 
major problem concerning them is that they have an infinite variance [3] and 
therefore their practical applicability is somewhat limited. Also they gen- 
erally do not possess a density expressible in an analytic form. This study 
proposes a dispersion measure for them, drawing ideas from Fisher informa- 
tion, differential geometry and most importantly, the uncertainty principle 
for Fourier transform pairs [5j. The study begins with a brief discussion on 
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characteristic functions and their relation to Fourier transforms and their 
properties, proceeds to a brief presentation of stable distributions and accu- 
mulates in defining a concept of characteristic curvature, which is proposed 
as a suitable measure of dispersion for class of stable distributions. 

2 Characteristic functions of random variables 
and Fourier transforms of probability mea- 
sures 

The following presentation will follow loosely [3]. In mathematics it is com- 
mon to obtain solutions to given problems by transforming the problem in 
some other space, solve the problem there, and then transform the solu- 
tion back to the original space. One of the most important transforms is the 
Fourier transform. It is an integral transform that transforms some function, 
or measure in so called frequency space. Fourier transform for a measure on 
R is defined as: 

Definition 1 Let F(x) be a measure on measurable space (R, £>r) ; where Br 
is the Borel a-algebra on R, then it has a Fourier transform defined as: 

<>oo 

F(u)= / e iux dF(x) (1) 

J — oo 

where e tux = cos(ux) + ism(ux). Moreover, if there exists a function f(x) 
such that dF(x)=f(x)dx, then it is equivalent to 

/oo 
e iu *f(x)dx (2) 
■oo 

Definition 2 Given probability distribution F(X < x), x G R, the charac- 
teristic function of the random variable X is defined as: 

/oo 
e iux dF{X < x) (3) 
■oo 

Since \e lux \ = 1, ipx{u) is defined for all u G R for every F(X < x). 

Definition Z If a probability distribution F(X < x) has a density f , its 
characteristic function is: 

/oo 
e iux f(x)dx (4) 
-oo 

That is, for probability distributions with density f(x), the characteristic 
function is just the Fourier transform of the density f(x). 
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We proceed with some properties of characteristic functions. 

Theorem 1 A characteristic function is uniformly continuous on R and sat- 
isfies the conditions 

(p x (0) = 1 ,\vx(u)\ < 1 (-00 <u<oo) (5) 

Proof: see JEI 

Theorem 2 Given random variables X and Y = aX + b, where a and b are 
constants, then the characteristic functions of the random variables X and 
Y are connected by the equation 

<py(u) = <p x (auy bu (6) 

and if a > 0, then their distribution functions satisfy the relation 

F(Y <x) = F(X < (7) 

a 

Proof. We note that 

ip Y (u) = E{e iuY ) = E{e iu{aX+b) ) = e iub E{e iuaX ) = e iub Vx {au) 
and for a > 

F(Y <x) = F(aX + b<x) = F(X < ^—^) 

a 

QED. 

Theorem 3 The characteristic function of a random variable X is real if 
and only if the distribution function of the random variable X is symmetrical, 
that is for every x the following equation holds: 

F(X <x)+ F(X <- x ) = l 

Proof. See [3]. The following theorem is of great importance: 

Theorem 4 (Convolution theorem) The characteristic function of the sum 
of two independent random variables is the product of the characteristic func- 
tions of the summands. That is, 

(Px+y{u) = ip x (u)(p Y (u) (8) 

Proof. Given X and Y independent, then e luX and e luY are also indepen- 
dent. Therefore, 

<p x+Y (u) = E(e iu{x+Y) ) = E(e iuX e iuY ) = E(e iuX )E(e iuY ) = ip x (u)ip Y (u) 

As expectation of the product of independent random variables is the product 
of the expectations, QED. 
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3 Stable probability distributions 

In this section, we introduce stable distributions and explore their most im- 
portant properties. We shall follow loosely [3]. 

Definition 4 The distribution function F(X < x) is called stable if to every 
a\ > 0, bi, a,2 > 0, b 2 there correspond constants a > and b such that the 
equation 

F(X < a x x + bi) * F(X < a 2 x + b 2 ) = F(X < ax + b) (9) 

holds. 

Equivalently, we can define stable distributions as: 

Definition 5 The distribution function F(X < x) is called stable if to every 
a\ > ; bi, ci2 > 0, &2 there correspond constants a > and b such that the 
equation 

F( < x) * F( <x) = F( < x) (10) 

a\ a 2 a 

holds. 
or, still: 

Definition 6 The distribution function F(X < x) is called stable if to every 
ai > 0, b\, a% > 0, b 2 there correspond constants a > and b such that the 
equation 

lf x-b (u) = (p x-b 2 (u)(fi x-b 1 (u) (11) 

holds. 

In words, a stable distribution F is such, that given random variable X 
and two affine transformations of it distributed with the same law F, then 
the sum of those transformed random variables is distributed according to 
some other affine transformation of X with the same law F. 

3.1 Canonical representation of stable laws 

Next theorem gives an explicit function class for the characteristic functions 
of stable distributions. 
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Theorem 5 In order that the distribution function F(X < x) is stable, it is 
necessary and sufficient that its characteristic function is: 

ipx{u) = e i«i (12) 

where a, 7, c > are constants (7 €R 7 — l</3<l,0<a<2 and 
u{u, a) = tan( z y) if a 7^ 1 and lj(u, a) — - In \u\ if a = 1 

The proof is long and can be found in [3]. Following is important for our 
purpose of this study: 

Corollary 1 A stable distribution function F(X < x) that is symmetric 
around zero^\, has the characteristic function 

<p x (u) = e~ c ^ a (13) 

This results from the properties of characteristic functions discussed above. 

3.2 Sums of symmetric, stable random variables 

Now we examine the summing of independent, identically distributed sym- 
metric, stable random variables. As the study of stochastic processes and 
therefore stochastic modelling is based on the properties of sums of random 
variables, in this chapter we study them in the case of symmetric, stable 
distributions. 

Consider the set of independent, identically distributed random variables 
Xi i G (1,2,..., t) with Xi having the characteristic function 

<p x .(u) = e- c ^ a \/i (14) 

Let us consider the sum of such random variables: 

t 

S t = J2 x i = x i+ x 2 + --- + x t (15) 

i=l 

We wish to determine the characteristic function of St- As the convolution 
theorem states that the characteristic function of the sum of iid. random 
variables is the product of respective characteristic functions, we shall have: 

t 

Vs t {u) = '[[<px i {u) (16) 
i=i 



this study, the expressions symmetric around zero and symmetric are considered 
synonymous 
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Substituting, 

tp St (u) = ( e - c M<7 = e -tc\n\ a ( 17 ) 
If we set C = C(t) = ct, we have: 

= (e~^y = e- CitMa (18) 

That is, the sum of t iid. random variables Xi has the same characteristic 
function as X* up to a linearly growing scaling factor C 

3.2.1 Example: The Gaussian case iV(0,a 2 ) 
Consider random variable Xi with characteristic function: 

m (u) = e-^ 2 (19) 

Now we have for the sum St the characteristic function: 

ip St {u) = e-W (20) 

a 2 (t) = aH (21) 
a(t) = o\Tt (22) 



If we set, 
or equivalently: 
we will have: 



cp St (u) = e -^W (23) 



As we can see, the variance of the sum grows linearly as a function of the 
number of summands. This kind of behaviour for sums of iid. random vari- 
ables, that is the linear growing of variance as a function of time (summands) 
is called normal diffusion. If one solves eg. the heat equation, one obtains a 
time-dependent normal distribution that has a linearly growing variance as 
a function of time [2]. 



3.3 Scaling of symmetric, stable distributions 

Next we concentrate on the scaling properties of stable distributions. As by 
definition the sum of t stable iid. random variables has the same distribution 
up to scale, the natural question is: what is that scale and how doest it 
evolve when the number of summands is increased! We answer this question 
straightforwardly by presenting the following theorem: 
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Theorem 6 (Scaling theorem) Given symmetric, stable random variable X 
with characteristic function <fx(u) = e~ c ' u ' Q and distribution function F(X < 
x). Then the sum S t oft iid. random variables X has the distribution: 

F(t°X < x) (24) 

or equivalently 

F(X < 4) (25) 

t a 

Proof. We have the characteristic function: 

= e~ c| " r (26) 

Suppose we scale the random variable X such that X — ► sX, where s > 0. 
The characteristic function of the scaled random variable sX will be: 



<Psx(u) = e~ c ^ a (27) 



Which is equivalent to: 



Vsx(u) = (28) 

In the previous section we showed, that the sum of t independent, identically 
distributed symmetric, stable random variables X has the characteristic func- 
tion: 

VSt {u) = e-^\ a (29) 

by comparing, we note that these two expressions are the same, if t = s a 
or s = t~ . Hence if we scale the random variable X by a factor t~ , the 
scaled random variable t^X has exactly the same distribution as the sum of 
t iid. random variables X. In distribution function notation, the sum of t 
independent, identically distributed symmetric, stable random variables X 
has the distribution: 

F{t^X < x) (30) 

or equivalently: 

F(X < J) (31) 

t" 

where F(X < x) is the distribution of X. QED. So the probability dis- 
tribution spreads out proportionally to i« when number of summands t is 
increased. This is the fractal property of sums of stable random variables, 
the sum distribution looks exactly the same as the original, it is just zoomed 
by a factor t~. This scaling property will come up later on in an interesting 
context. 
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3.4 Some remarks on the scaling properties 

As we have seen, stable distributions obey a scaling law of the type t« for 
sums of t iid. symmetric, stable random variables. We have also seen that 
the scale parameter c > in the characteristic function grows linearly as a 
function of the number of summands (or time, if you will). These properties 
do tell us a whole lot about the sum behaviour of stable random variables. 
In the Gaussian case one can directly state how the standard deviation grows 
over time if we consider the sum as a 1-dimensional random walk. This is 
extremely useful, as for example in the context of mathematical finance, stan- 
dard deviation of the underlying price process is called volatility. However, 
as the Gaussian case is the only case where variance is finite, we can not give 
predictions of the growth of standard deviation over time in other cases when 
a < 2. Basically we face the problem that we do know the time-evolution of 
the dispersion of the sum-distribution, but we do not know how to measure 
that dispersion in general. This is the topic of the next chapter. 

4 Characteristic curvature 
4.1 Motivation 

We next proceed to propose a measure for dispersion for symmetric around 
zero stable distributions with a G [1,2]. The idea is to think of a measure 
that somehow characterises the dispersive properties of the characteristic 
functions. We know that characteristic functions determine probability dis- 
tributions completely [3], whether the moments are finite or not. So the 
general dispersive properties of different probability distributions must be 
somehow incorporated in the properties of characteristic functions. We know 
from Fourier analysis, that the more spread the original function, the more 
concentrated its Fourier transform and vice versa [5]. Because characteristic 
function is the Fourier transform of the probability measure, we can apply 
this idea here. It is therefore natural to think that if the probability dis- 
tribution is greatly dispersed in space, the characteristic function should be 
condensed around the origin in the frequency space. This concentration of 
characteristic functions around the origin is intuitively analogous to that the 
second derivative in its absolute value should be relatively large around the 
origin, because the characteristic function satisfies y?x(0) = 1 and it should 
decay fast when moved from the origin into either direction. Hower, be- 
cause the characteristic function here is symmetric around the origin, we can 
constrain to analyse the positive part of u. Therefore the negative of the 
second derivative should be relatively large for small u. To measure this con- 
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centration analytically, let us consider the weighted average of the negative 
second derivatives of ln-transform of (px{u) (assuming ipx( u ) monotonically 
decreasing Vw > 0) in the continuum limit, weights given by the charac- 
teristic function. That is, the integral transform that we shall call as the 

characteristic curvature: 

Definition 7 The characteristic curvature for symmetric around zero prob- 
ability distribution F(X < x) is 

V27T Jo ou l 

where (px{u) is the characteristic function of the distribution and u > 0. 

The constant -j= is suitable for fixing the proper scale, as we shall see later 
on. Note that this integral resembles the definition of Fisher information, 
differences being that now the likelihood function is the characteristic func- 
tion, and the parameter is u. By performing the derivation and simplifying, 
the characteristic curvature can be formulated equivalently as: 

$ = ^L r (^j^ - V ' x (u))du (33) 
V27T Jo <Px{u) 

which resembles some curvature concepts in differential geometry. For the 
sake of comparison, the standard deviation of a probability distribution with 
characteristic function (px(u) is given as: 



a = yJ-WxM) 2 ~ MM (34) 

see eg. [3]. 



4.2 Characteristic curvature for symmetric, stable dis- 
tributions 

Theorem 7 Given symmetric around zero stable distribution with charac- 
teristic function (fx{u) = e~ c ' M ' a its characteristic curvature satisfies 

*(c; „) = * r { - ffh /5 W K^«fa = -H=aT(2-I)ci = ^r(2 
V2^Jo du 2 v a 1 H\F2sk v 

(35) 

where T(2 — ^) is the Euler Gamma function and H is the Hurst exponent. 
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Proof. 

Let us begin with the characteristic function of a symmetric stable distribu- 
tion for random variable X. 

Vx ( u ) = e - c l"l" uGK,c>0,0<a<2 (36) 

Because the characteristic function is symmetric around zero, we consider 
the positive half-space representation of it: 

<p x ( u ) = e' cua u>0, c> 0, < a < 2 (37) 

Taking logarithms, we will get 

\nip x (u) = -cu a (38) 

Derivating twice with respect to u: 

= -a(a - l)cu a -' 2 (39) 



d 2 In ip x (u) _ , -2 



du 2 

Which is the same as 

d 2 In lpx(u) , x „ 9 , . 

— = a{a ~ l)cu (40) 

Substituting this into characteristic curvature, we shall have: 

$(c; a) = -= / {a(a - l)cu a - 2 }e' cua du (41) 

V27T Jo 

Which is equivalent to: 

2 f°° 
$(c; a) = ~7=oi{a - l)c J u a - 2 e~ cua du (42) 



Next we calculate the integral: 

PCX) 

/ u a - 2 e~ cua du (43) 
Jo 



Substituting t = cu a dt = cau a 1 du du = ca ^ a -i The integral will 
become: 



which is 

roc 



/•oo a-2 -t 

/ r rft 44 

7 caw"- 1 

u _1 e"*dt (45) 



1 
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Noting that t = cu a 44> u = ^- the integral becomes 



ta 

i r t- 



ca 



-e- l dt (46) 



C <* 



Taking the constant outside the integral gives: 

1 f°° 

— / r^e-'dt (47) 

acc a J 

By simplifying and adding zero to the exponent of t, we shall have: 

i-r / t( 1 -=)- 1 e-*dt (48) 



ac a Jo 

Recalling that the Euler Gamma function is defined as 

POO 

r(z) = / t z - l e- l dt zeC (49) 

Jo 

We will obtain the elegant result: 

u «-2 e -cu" du = 2llr(l - -) (50) 
a a 

Substituting this to the characteristic curvature, we will obtain: 

2 c~~ l 1 

&cr,a) = —=a(a-l)c——T(l--) (51) 
V27T ol a 

Simplifying, it will become: 

$( c;a ) = (a _ 1)c i r(1 _ 1 } 
V27T a 



Because the Gamma function satisfies 

Y(z + l) = zY(z) (53) 

We have 

r ^ 2 -b = ^~ b^ 1 - b * r(1 ~b = 7^ T{2 ~ b ^ 

a a a a a — I a 

Substituting this in the characteristic curvature: 

$(c; a) = -^=(a - l)c« — —T(2 --) = ^=ac«T(2 - -) (55) 
V27T ol — 1 a V27T a 



We will finally get the most compact form for characteristic curvature: 



2 1 

$(c; a) = ^=aV(2 - -)c- (56) 
V27T a 



QED 



4.3 Properties of characteristic curvature for stable distributions 
and some restrictions 12 



4.3 Properties of characteristic curvature for stable 
distributions and some restrictions 

We have shown that the characteristic curvature for stable distributions sat- 
isfies: 

2 1 i 

$(c; a) = ^=«r(2 )c» 

V27T a 



As the gamma function has singularities at x — 0, —1, —2, —3..., characteris- 
tic curvature diverges when a approaches \. If we restrict that 1 < a < 2, 
we will have following desirable properties for characteristic curvature: 

Theorem 8 Given 1 < a < 2 and Vc > ; the characteristic curvature 
$(oj, c) satisfies the following properties: 



(1) $(c;a)>0 

(2) > 

Proof. (1) is trivial. We prove hence only (2). By derivating: 



(9$ 2 , 1, l-a 

— = -=r(2 - -) c — > o 

QED. In other words, characteristic curvature is strictly increasing in c, which 
is the scale factor, and is always positive. 

4.4 Characteristic curvature for Gaussian distribution 

Let us calculate the characteristic curvature for centered Normal distribution 
with variance a 2 . The characteristic function is 

<p x {u) = e-^ 2u2 

Equivalently a = 2 and c=\o 2 . By substituting these values, we will get: 



Because T(|) = we will obtain 



1 2, 4/1 2V ^F 2^ 1 
$(2, -a ) = -^=\ -a 2 -^- = —j=\ -a 2 

1 '2 ; 2 2 2 



Which simplifies to: 



$(2, \a 2 ) = V2^j\a = a 



4.5 Characteristic curvature for Cauchy distribution 
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This is very fascinating result! The characteristic curvature of a centered 
Gaussian with variance a 2 is just the standard deviation a. Note that also 
holds: 

$ 2 (2,V)=a 2 (57) 

, which is obvious, of course. 



4.5 Characteristic curvature for Cauchy distribution 

Let us calculate the characteristic curvature for Cauchy distribution, that is 
a = 1 and c is unspecified. The characteristic function is: 

Mu) = e" cH 
Substituting these values we will have: 

•d.c) = ^r(i) 



Because T(l) = 1, we will have simply: 

2c 



$(l,c) 



2tt 



4.6 Sums of symmetric, stable random variables and 
characteristic curvature 

We proceed to consider the characteristic curvature for sums of independent, 
identically distributed symmetric, stable random variables Xi with charac- 
teristic function (fx,( u ) = e~'"'° Vi G (1,2,3, ...,t). We recall that the sum 

t 

S t = ^ = X l + X 2 + ■ ■ ■ + X t (58) 
i=i 

has the characteristic function 

¥>*(«) = e-*M a (59) 
We calculate the respective characteristic curvature: 

2 1 i 

$(tc; a) = ^=oT(2 )(tc)° (60) 

V27T a 



4.6 Sums of symmetric, stable random variables and 
characteristic curvature 
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which becomes 

2 lii 

$(tc; a) = -=aF(2 )c°t° (61) 

V27T a 

which equals 

$(tc; a) = $(c; a)t * (62) 



This is very elegant result! The characteristic curvature scales like ~ 
which is exactly the same scaling behaviour that stable random variables 
possess in general. 

4.6.1 Example: X t ~ N(0, a 2 ) 

Given that Xj ~ iV(0,<7 2 ) Vi That is, we have c = |a 2 , a = 2, so for the 
sum of t such random variables, we have: 

$(tia 2 ;2) = $(ia 2 ;2)^ (63) 

We recall from earlier results that $(|cr 2 ; 2) = a, therefore we obtain: 

<5>(t-(T 2 ; 2) = aVi (64) 

Which is exactly the same result we had earlier for normal diffusive be- 
haviour. Note that the square of the characteristic curvature is additive in 
the Gaussian case , that is when it is the variance, as: 

$ 2 (tc; a) = $ 2 (c; a)t« = crH (65) 

4.6.2 Example: X { ~ Cauchy (c) 

Let Xi have a Cauchy distribution with some scale factor c. We have a — 1, 
so for the sum of t such random variables, we have: 

$(te; 1) = $(c; 1)*t (66) 

recalling earlier results, we then have: 

$(tc; 1) = ^=t (67) 



which is the defining characteristics of so called ballistic diffusion. 



5 Conclusions 



15 



5 Conclusions 

We have proposed a dispersion measure for stable distributions with 1 < 
a < 2. Characteristic curvature has very intuitively appealing roots and 
interesting connections to various branches of mathematics. The scaling 
property; the familiar results of superdiffusive characteristics and the reduc- 
tion to standard deviation in the Gaussian case are interesting implications 
of the definition. The abundance of Gamma function might lead to other in- 
teresting connections of stable distributions and e.g. Riemann zeta function 
as they are intimately related to each other. Obviously, more research is vi- 
tally needed to establish general results for the applicability of characteristic 
curvature for general probability distributions. Also, one interesting field of 
applied research would be that of option pricing for stable price processes. 

References 

[1] A. Ya. Khintcine, Paul Levy Sur les lois stables, C. R. Acad. Sci. Paris 
202, 374-376 (1936) 

[2] L. C. Evans Partial differential equations, American mathematical soci- 
ety, ISBN 0-8218-0772-2 (1998) 

[3] B. V. Gnedenko, A. N. Kolmogorov Limit distributions for sums of in- 
dependent random variables, Addis on- Wesley (1954) 

[4] J. Jacod, P. Protter Probability essentials, Springer ISBN 3-540-43871-8 

(2004) 

[5] H. Weyl Theory of groups and quantum mechanics, Dover, NY (1950) 

[6] B. B. Mandelbrot Fractals and scaling in finance, Springer ISBN 0-387- 
98363-5 (1997) 



